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We introduce and analyze a new quantity, the path integral ideal, governing the flow of generic 
discrete theories to the continuum limit and greatly increasing their convergence. The said flow 
is classified according to the degree of divergence of the potential at spatial infinity. Studying the 
asymptotic behavior of path integral ideals we isolate the dominant terms in the effective potential 
that determine the behavior of a generic theory for large discrete time steps. 
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Path integrals present a rich and flexible formalism for 
dealing with quantum and statistical theories Q, Q that 
has proven extremely useful for handling symmetries, de- 
riving non-perturbative results, establishing connections 
between different theories Q, and extending the quan- 
tization procedure to ever more complicated systems. 
They have served as catalysts for the exchange of key 
ideas between different areas of physics, most notably 
high energy and condensed matter physics Q. Today, 
analytical and numerical approaches to path integrals 
0, Q play important roles not only in physics but also 
in chemistry and materials science, and are acquiring a 
prominent role in mathematics and modern finance |8j . 

Further development of the path integral method is 
constrained by the small number of solvable models, as 
well as by our rather limited knowledge of their precise 
mathematical properties. In fact, most of our knowledge 
is negative, e.g. we know which trajectories do not con- 
tribute to the path integral rather than which do. One of 
the few positive statements concerning path integrals is 
that relevant trajectories exhibit stochastic self-similarity 
PJ. As a result they have non-trivial fractal dimension 
and jaggedness [{| . Researchers working on numerical ap- 
proaches to path integrals have successfully utilized these 
kinematic consequences of self-similarity to produce effi- 
cient path-generating algorithms 0. 

In a recent series of papers we have investigated 

the dynamical implications of stochastic self-similarity by 
studying the relation between discretizations of path in- 
tegrals with different coarseness. This has resulted in a 
systematic analytical procedure that may be used to re- 
duce path integral error to 0(e p ) for arbitrary p 6 N, 
where e is the discrete time step. Note that e = T/N, 
T being the time of propagation and N the discretiza- 
tion coarseness. This reduction of error brings about a 
substantial increase in the speed of numerical algorithms. 
Self-similarity played a crucial role in this procedure in 
that it allowed us to derive an integral equation relating 
discretized theories viewed at different coarseness and to 
solve it in terms of an asymptotic series. The asymptotic 
expansion, however, implies that the obtained method is 
directly applicable only for e < 1. 

The fact that we can arbitrarily decrease the error 
points to the possibility that one can extend the for- 



malism and obtain exact information (i.e not given as 
a power series in e, and so valid even for large values 
of e) about the continuum theory. Large e corresponds 
to long times of propagation, precisely what interests us 
in quantum field theory (or in modern finance). Equiv- 
alently, in condensed matter and materials science this 
corresponds to the physically most interesting region of 
small temperatures. Large e behavior is also central for 
determining the energy spectrum of a given model, and 
as such is applicable in many areas of physics y| (e.g. 
atomic and molecular physics, quantum dots). 

In this paper we introduce and analyze a key new quan- 
tity, the path integral ideal, which governs the flow of 
generic discrete theories to the continuum. We show that 
this flow is classified according to the degree of diver- 
gence of the potential at spatial infinity. We isolate the 
dominant part of the effective potential and show that 
it mediates the large e behavior of the theory. In this 
way we extend the previous method and derive the first 
properties valid for large e. 

In we showed that at the core of the new method 
one finds a generalization of Euler's summation formula 
to path integrals. The Euler sum formula for the or- 
dinary integral fT dtf(t) can be viewed as an explicit 
construction of an ideal discretization f*(t;e) of the 
integrand, where e = T/N is the discrete time step, 
and /*(t;e) — * f(t) in the continuum limit e — * 0. 
The discretization is ideal because for any coarseness N 
it exactly reproduces the integral we seek to evaluate 

E^=i £ /*(*n;e) = Sq dtf{t), where t n = ne. The Euler 
sum formula gives f*(t; e) as an infinite power series in e. 
The result is /* = /-§/ - ^/+ • ■ ■■ Truncation of /* 
to the first p terms depends on f(t) and its first (p — 1) 
derivatives and determines the integral up to 0(e p ). 

The generalization of this to path integrals showed 
the existence of an ideal discretized effective action S* 
having the property that at any discretization coarse- 
ness N its discretized amplitudes exactly reproduce the 
corresponding continuum amplitudes calculated with the 
starting action S. For (Euclidean) actions of the form 
S = J Q dt (| q 2 + V(q)), the ideal discretized effective 
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action was shown to be 

S *=Y, (f e +eW*(S n ,q n ;e)) , (1) 

n—o ^ ' 

where S n = q n +i — q n are the discretized velocities and 
I 



q n = i + q n ). For obvious reasons we call the func- 
tion W* (5, q; e) the path integral ideal. We have shown 
[lll | that this function is even in the first argument, that 
in the continuum limit it goes over to the potential V(q), 
and that it satisfies a novel integral equation relating dis- 
cretizations of different coarseness: 



r, p + ao 
7T£ ./ o 



exp 



2 V 2 



o ye 
M+ 4 + 2 ; 2 



In previous papers we asymptotically expanded the in- 
tegral in the above equation obtaining a solution in the 
form of an infinite power series in e, namely 



W*{5,q;e) 



n=0 
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(3) 



The coefficient functions C n (S, q) were found to be given 
in terms of algebraic equations. Using MATHEMATICA 
5.0 we obtained explicit solutions for the first nine coeffi- 
cients. The truncation of W* to the first p terms depends 
on the potential V and its first (2p — 2) derivatives and 
allows us to calculate the path integral up to an 0(e p ) er- 
ror term. This has brought about a substantial speedup 
in the numerical calculations of path integrals. 

Our aim is to study the asymptotic properties of the 
above sum. Similar asymptotic expansions often arise in 
physics. For example, perturbative expansion of path in- 
tegrals (expansion in some small coupling constant g) 
generally leads to infinite sums of the form J2 n A n g n 
where, for large n, the coefficients go as A n ~ (— ) n n\. 
Therefore, in analyzing the asymptotic properties of per- 
turbative expansions one typically encounters sums like 
Yl n (~ ) n n\g n . Using Stirling's formula we find that the 
n-th term in the sum is asymptotically (— gn) n . Hence, 
the sum is divergent since its terms start increasing af- 
ter n c ~ l/g. The divergence is a result of a cavalier 
interchange of infinite sums and integrals at the very 
beginning of the asymptotic expansion procedure. Al- 
though divergent, asymptotic sums give excellent approx- 
imations (for small g). The best approximations are ob- 
tained by truncating the sum at n c . In principle one can 
often do much better than that. For example, leading 
strong coupling behavior may be obtained using Borcl 
resummation jj]. More sophisticated techniques based 
on the variational perturbation method allow us even to 
extract sub-leading strong coupling behavior p| . This is 
the ultimate reason why we are interested in the asymp- 
totic properties of path integral ideals, since, as it stands, 
the sum in eq. © can only be used for small e. Exact, 
in this sense, means valid for all e. 

It is instructive to begin our investigation of the 
asymptotic properties of ideals by looking at the cor- 



responding expressions at the p = 3 level. In fact, 
rather than look at the whole ideal we first focus on 
its velocity independent part - the effective potential 
V*(q; e) = W*(0, q; e). At the p = 3 level we found [13 



V* = V 




(4) 



Note that, as a consequence of the V term being neg- 
ative, it is possible for the above effective potential to 
be unbounded from bellow. The top left plot in Fig. 1 
shows the p = 3 level effective potential for the case of an 
oscillator with quartic coupling V\ (q) = \q 2 + \ x g q 4 - We 
see that the corresponding effective potential has critical 
points ±g c , where q c ~ Problems arise when using 
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FIG. 1: (Color online) Effective potentials at level p = 3 (solid 
thick line) and their critical points for the case of four different 
models: oscillator with quartic coupling Vi(q) = + q 4 
(the case e = 0.5, g = 1 is shown top left); particle moving in a 
Morse potential V 2 {q) = A(e~ 2Xq - 2e~ Xq ) (the case e = 0.1, 
A = 10, A = 2 is shown top right); particle moving in a 

modified Poschl- Teller potential V 3 (q) = \ fj^' 1 ^ (the case 
e — 1, a — 0.5, j3 = 1.5 is shown bottom left); particle in a 
Lennard- Jones potential T4(r) = — ^ (the case e = 0.01, 
C = 1, D = 20 is shown bottom right). For comparison the 
original potentials have been plotted with a dashed thin line. 
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this effective action in the path integral to integrate over 
trajectories with \q\ > q c . The top right plot in Fig. 1 
shows the effective potential for a particle moving in a 
Morse potential V 2 {q) = A(e~ 2Xq - 2e~ Xq ). In this case 
there is only one critical point q c ~ These two ex- 
amples illustrate a general property: effective potentials 
become unbounded from bellow at large positive (nega- 
tive) values of q if the original potential diverges at posi- 
tive (negative) spatial infinity. For example, the modified 

does not di- 
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Poschl- Teller potential Va(q) 

verge at infinity and leads to a well behaved effective 
potential (bottom left plot in Fig. 1). Another interest- 
ing class of potentials are those that diverge at a point. 
From eq. Q we see that potentials that diverge at a point 
go as (q — qo)~ n have no critical points for n < 2, the 
best example of this being the Coulomb potential. The 
Lennard- Jones potential Vi(r) — C/r 12 — D/r 6 diverges 
faster at the origin giving its p — 3 effective potential a 
critical point (bottom right plot of Fig. 1). 

Unboundedness from bellow is specific to certain p lev- 
els only, e.g. p = 5 level effective potentials are well be- 
haved for all models. At first hand it is surprising to 
have unboundedness from bellow for any model at any 
p level. To make sense of these effective potentials one 
needs to cut off coordinates at the critical points. In fact, 
far from surprising, this is just the analogue of truncating 
asymptotic series to a critical number of terms n c . 

In solving an ordinary integral in an asymptotic ex- 
pansion about a small parameter g one obtains an infi- 
nite sum. In many cases this sum is not convergent but 
asymptotic and the best approximation is obtained by 
truncating it at n c . As already shown n c typically be- 
haves as 1/g, or as some power of this, depending on the 
precise integral being evaluated. In deriving the integral 
equation for ideals we asymptotically expanded path in- 
tegrals (the discrete time step e playing the role of small 
parameter) and obtained the result in the form of an ordi- 
nary integral. We have now come to the conclusion that 
in many cases this integral is not convergent but asymp- 
totic (depending on the behavior of the original potential 
at infinity) and that the best approximation is obtained 
when the effective potential (and ideal) are truncated at 
q c . Typically q c behaves as 1/e, or some power of this 
(depending on the model we are working with) . 

Let us point out that the original Monte Carlo sim- 
ulations used to check the analytically derived results 
|10L 1 1 1| worked well without any need for truncation at 
q c . The reason is that the path generating procedures in 
numerical simulations effectively create their own cut- 
offs. The width of generated trajectories is typically 
a = e/2. The requirement that trajectories do not 
overshoot the critical point is that q c is much larger than 
a. For the oscillator with quartic coupling this is achieved 
for N > (g/2)^ 3 T. For the case g = 10 and T = 1, for 
example, this criterion is satisfied even for the smallest 
A's (roughest discretizations) and in those simulations 
the critical points were never reached. 

Asymptotic behavior of ideals is determined by the 



places where the original potential diverges as well as the 
degree of that divergence. Finite potentials lead to well 
behaved effective potentials and ideals. An illustration 
of this is the modified Poschl- Teller potential. Similarly, 
the Morse and Lennard- Jones potentials have path in- 
tegral ideals that are well behaved at right infinity (i.e. 
where the original potentials do not diverge). 

In the remainder we will concentrate on potentials 
whose only divergence is at spatial infinity. If the domi- 
nant behavior of V at infinity is q n then from eq. we 
see that this is also the dominant behavior of the posi- 
tive part of the p — 3 level effective potenital, while the 
negative part goes as e 2 q 2n ~ 2 . The effective potential is 
well behaved for n < 2. The case n = 2 is integrable 
and so need not interest us here. Still, let us note that 
it too is well behaved provided e is small enough. The 
case n > 2 produces a critical point. The simplest way to 
estimate the critical point is to note that in its vicinity 
the positive and negative terms in the effective potential 
balance out. Therefore, for potentials that diverge with 
a power n > 2, we have q c ~ l/e 2 ^ n ~ 2 \ Similarly, for 
potentials diverging exponentially we find q c ~ — lne 2 . 

From studying the effective potential at the p = 3 level 
we have learned which potentials give well behaved effec- 
tive actions (and ideals) and which have critical points, 
and have evaluated these critical points from the way the 
original potentials diverge at infinity. We have also seen 
that it is important to study the ideal in the vicinity of 
the critical points, where the great majority of its terms 
are subdominant. It is not difficult to see that the domi- 
nant part of the complete p = 9 level solution [10] equals 
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We have written the dominant part of the ideal as V£ om 
as it is velocity independent and part of the effective po- 
tential. We have numerically checked that for all models 
considered this gives an excellent approximation to the 
ideal even far from critical points. From the above for- 
mula we see that effective potentials are well behaved 
(bounded from bellow) for p = 1,2,5,6,9,10,..., it is 
only at the remaining p levels that one may find critical 
points. This is the first indication that critical points are 
artifacts of the asymptotic expansion procedure and are 
not present in the exact path integral ideal. 

The terms in the above expression can be determined 
directly (up to numerical factors) without first solving for 
the ideal. To do this we introduce the following index: 
index(a&) = index(a) + index(6), mdex(V^) = 2 + a, 
index(e) = index(<5 2 ) = —2. Note also that the difference 
W* — V depends only on the derivatives of V. It follows 
that index(VF*) = 2, i.e all the terms in the ideal have 
the same index, as can be shown directly from eq. ©• 
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To determine the dominant terms at spatial infinity we 
need to: (a) look at terms that are products of as many 
V's as possible, (b) have as few derivatives as possible, (c) 
pick out terms with smallest power of e. All of this needs 
to be done with the constraint that the index remains 
fixed and equal to 2. The fact that the dominant part of 
the ideal is velocity independent directly follows from the 
relation mdex(S 2 ) = — index(V). Every added power of 
5 2 (ideals are even in the velocities) removes one power of 
V from the term and so decreases its order of divergence. 
One can also show that leading terms proportional to 
£ 2m+i are su bdominant to leading e 2m terms. As a result, 
the dominant part of the ideal is a power series in e 2 
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Leaving out for the moment the proof that the signs in 
the above series alternate we focus on the asymptotic be- 
havior of the coefficients c k which consist of products of 
(k + 1) V's with 2k derivatives. If V diverges as q n then 
c k ~ A k g( fe + 1 ) n— 2fc , where A k are constants that can't be 
determined from the index alone. From eq. |J5J it follows 
that the A k s are smaller than one and decrease with k. 
At q c we find c k e 2k ~ A k e -2n/("-2) j i. e . all the k de- 
pendance is in the constants A k . A similar result follows 
when V ~ e q . Therefore at q c the ideal is proportional to 
J2k(~ ) h A k - This is finite since it is smaller than ^ k {—) h 
which is regularized to 1/2. 

The toy model for the asymptotics of ideals (when all 
the Ak's are equal to 1) can be solved exactly. For po- 
tentials that diverge as q n the path integral ideal at a 
point Xq c equals W t * oy ~ \™ e -^/(n-2) -£ k (-f A^" 2 ^ = 
X n e -2n/{n-2) j(i + \n-2y We nnd that the alternating 
sign in eq. JSJ plays a crucial role in making this conver- 
gent at the critical point (A = 1). From W* < W^ oy it 
follows that, unlike its truncations at a given p level, the 
exact path integral ideal remains well behaved even at 
the critical points. We conclude that critical points, un- 
boundedness from bellow and truncations of coordinates 
are all artifacts of the asymptotic expansion procedure 



and are not present in the exact ideals. 

The only thing that remains is to prove that Vf om is 
indeed given in terms of an alternating sum. To do this 
we first set S = in eq. J2J , and use the fact that the ideal 
is even in the first argument. Next, we set q to be very 
large, a regime where the ideal may be approximated 
by its dominant term. In this regime q 3> y ~ y/e so 
that Vj om (q + y/2; e/2) is well approximated by the first 
two terms in its Taylor expansion about q. Finally, after 
calculating the remaining Gaussian integral we find 
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From this we can directly derive the first two terms of 
eq. (J5J. To get all the V d * om terms correctly it would be 
necessary to solve its integral equation asymptotically. 
This is a much simpler job than was solving eq. f2"|. Still, 
eq. (JJJ is good enough for us to extract three important 
properties of V d * om : (a) it is given as a power series in e 2 , 
(b) the terms in this series alternate in sign, (c) and the 
difference V£ om — V is proportional to V . 

To conclude, we studied the asymptotic properties of 
path integral ideals - expressions that allow us to con- 
struct the best discretized effective action for a general 
quantum theory. Path integral ideals govern convergence 
flow to the continuum limit allowing us to: (a) system- 
atically accelerate numerical algorithms, (b) investigate 
the behavior of quantum theories for large times of prop- 
agation (low temperatures) using extremely coarse dis- 
cretizations. We have shown that ideals are classified ac- 
cording to the degree of divergence of the original poten- 
tial at spatial infinity, and have isolated dominant terms 
responsible for the behavior of the ideals of a generic 
theory for large discrete time steps. We have shown that 
exact path integral ideals are always well behaved, in 
contrast to the behavior of some of their e p truncations. 
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